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Abstract 

Flameless  combustion  offers  many  advantages  over  conventional  combustion,  particularly  uniform  tem¬ 
perature  distribution  and  lower  emissions.  In  this  paper,  a  new  strategy  is  proposed  and  adopted  to  scale  up 
a  burner  operating  in  flameless  combustion  mode  from  a  heat  release  density  of  5.4-21  MW/m3  (thermal 
input  21.5-84.7  kW)  with  kerosene  fuel.  A  swirl  flow  based  configuration  was  adopted  for  air  injection 
and  pressure  swirl  type  nozzle  with  an  SMD  35-37  pm  was  used  to  inject  the  fuel.  Initially,  flameless  com¬ 
bustion  was  stabilized  for  a  thermal  input  of  21.5  kW  (Q'"=531  MW/m3).  Attempts  were  made  to  scale  this 
combustor  to  higher  intensities  i.e.  10.2,  16.3  and  21.1  MW/m3.  However,  an  increase  in  fuel  flow  rate  led  to 
incomplete  combustion  and  accumulation  of  unburned  fuel  in  the  combustor.  Two  major  difficulties  were 
identified  as  possible  reasons  for  unsustainable  flameless  combustion  at  the  higher  intensities,  (i)  A  constant 
spray  cone  angle  and  SMD  increases  the  droplet  number  density,  (ii)  Reactants  dilution  ratio  (Rdu) 
decreased  with  increased  thermal  input.  To  solve  these  issues,  a  modified  combustor  configuration,  aided 
by  numerical  computations  was  adopted,  providing  a  chamfer  near  the  outlet  to  increase  the  Rdij.  Detailed 
experimental  investigations  showed  that  flameless  combustion  mode  was  achieved  at  high  intensities  with  an 
evenly  distributed  reaction  zone  and  temperature  in  the  combustor  at  all  heat  intensities.  The  emissions  of 
CO,  NOv  and  HC  for  all  heat  intensities  (0  =  1-0.6)  varied  between  1 1-41,  6-19  and  0-9  ppm,  respectively. 
These  emissions  are  well  within  the  range  of  emissions  from  other  flameless  combustion  systems  reported  in 
the  literature.  The  acoustic  emission  levels  were  also  observed  to  be  reduced  by  8-9  dB  at  all  conditions. 
©2014  The  Combustion  Institute.  Published  by  Elsevier  Inc.  All  rights  reserved. 
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1.  Introduction 

Flameless/Mild  combustion  has  gained  signifi¬ 
cant  importance  due  to  its  ability  to  suppress  ther¬ 
mal  NO  formation  and  improve  thermal  efficiency 
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of  combustion  systems.  Flameless  combustion  has 
been  primarily  identified  with  gaseous  fuels  and 
extensive  work  has  been  reported  1—10].  Scaling 
the  flameless  combustors  to  higher  intensities  has 
been  proposed  in  recent  studies  reported  in  the  lit¬ 
erature  [7,8,11].  A  brief  summary  of  various  high 
intensity  flameless  combustion  systems  with  gas¬ 
eous  fuels  is  listed  in  "able  1.  Liickerath  et  al. 
[11  have  developed  a  Forward  Flow  (FF)  com¬ 
bustor  configuration  with  a  thermal  input  of 
475  kW  and  heat  intensity  of  240  MW/m3  (at 
20  bar).  Kumar  et  al.  have  scaled  up  a  high- 
intensity  combustor  (5-150  kW  thermal  input) 
with  new  scaling  methodology  and  compared  var¬ 
ious  existing  scaling  techniques  i.e.  Constant 
Velocity  (CV),  Constant  Residence  Time  (CRT) 
and  Cole’s  approach  with  the  proposed  technique. 
The  comparison  of  Weber  12]  shows  that  CRT 
approach  is  relatively  better  for  scaling  swirl  type 
combustor  configurations.  They  have  hinted  at 
the  need  of  maintaining  high  reactant  dilution 
rates  to  ensure  that  flameless  combustion  mode 
is  achieved  in  scaled  combustors.  These  types  of 
combustor  configurations  with  high  heat  intensity 
are  expected  to  be  useful  in  gas-turbine  applica¬ 
tions.  Arghode  and  Gupta  [8]  have  demonstrated 
a  laboratory  scale  combustor  achieving  colorless 
distributed  combustion  with  a  high  intensity  of 
453  MW/m3  ( Qth  =  6.25  kW)  with  a  combustor 
volume  of -13  cm  and  Reverse  Flow  (RF)  con¬ 
figuration.  Scaling  of  these  concepts  with  low 
thermal  input  and  high  heat  intensity  render  the 
systems  very  complex  and  making  their  implemen¬ 
tation  highly  challenging.  Further,  very  little 
literature  is  available  in  the  field  of  scaling  of 
flameless  combustors  with  liquid  fuels.  Some  basic 
studies  on  flameless/mild  combustion  with  liquid 
fuels  [3,5,13-15  have  been  reported  recently. 
Traditional  industrial  burners  and  stationary 
gas-turbine  combustors  operate  with  liquid 
fuels  at  higher  thermal  inputs  (—1  MW)  and 
higher  heat  intensities  (100  MW/m3).  Therefore, 


Table  1 

Variation  of  heat  intensities  reported  in  literature  (SJ: 
Straight  Jet,  Forward  Flow,  FF  (i.e.  reactants  enters 
from  one  side  and  products  leave  from  opposite  side), 
Reverse  Flow,  RF  (i.e.  reactants  and  products  from 
same  side  of  the  combustor),  Qth\  Thermal  input  (kW) 
and  Q"'\  Heat  intensity  (MW/m3),  S:  Solid,  L:  Liquid,  G: 
Gas. 


Refs. 

Qth 

0" 

Fuel 

Remarks 

[2] 

5 

10 

G 

SJ,  FF 

[3] 

580 

0.024 

G,L,S 

SJ,  RF 

[5] 

0.4 

0.58 

G,  L 

SJ,  RF 

[6] 

20 

0.44 

G 

SJ,  RF 

[7] 

150 

5.6 

G 

SJ,  FF 

[8] 

6.25 

453 

G 

SJ,  RF 

[9] 

15 

0.3 

G 

SJ,  RF 

[11] 

474 

240 

G 

SJ,  FF 

additional  studies  are  required  to  investigate  the 
issues  related  to  the  scaling  of  high  intensity 
flameless  combustors  with  liquid  fuels  and  their 
relation  with  spray  characteristics. 

In  this  study,  a  swirl  based  combustor  operat¬ 
ing  in  flameless  combustion  mode  with  kerosene 
at  21.5  kW  (base  case,  0"'=5.37  MW/m3)  [16],  is 
developed  and  scaled  up  to  operate  at  85  kW 
(21.1  MW/m3).  Attempts  were  made  to  achieve 
flameless  combustion  with  higher  intensity  of 
10.2,  16.3  and  21.1  MW/m3  using  the  base  case 
combustor  configuration  with  increased  thermal 
input  of  40.8,  65.1  and  84.7  kW  respectively. 
However,  the  existing  combustor  configuration 
was  unable  to  achieve  stable  flameless  combus¬ 
tion.  Computational  and  experimental  studies 
were  carried  out  to  identify  the  causes  preventing 
successful  scale-up  of  the  combustor.  The  com¬ 
bustor  configuration  was  modified  by  providing 
a  chamfer,  thus  increasing  the  degree  of  recircula¬ 
tion  in  the  combustor,  allowing  stable  flameless 
combustion.  Computational  studies  show  that 
with  increased  chamfer  radius  ( Rc ),  recirculation 
of  the  combustion  products  and  fuel  residence 
time  increased.  Three  different  Rc  values  were 
considered,  for  the  10.2,  16.3  and  21.1  MW/m3 
cases  respectively,  and  shown  to  achieve 
combustion  with  low  emissions.  The  influence  of 
spray  characteristics  on  scaling  the  combustor 
was  studied  and  the  results  are  presented  in  this 
paper. 

2.  Computational  studies 

2.1.  Geometry  design  methodology 

The  base  combustor  was  designed  to  stabilize 
high  intensity  flameless  combustion  using  conven¬ 
tional  liquid  fuels.  Stabilization  of  flameless 
combustion  with  liquid  fuels  depends  on  three 
important  parameters, 

1.  Sauter  Mean  Diameter  (SMD)  of  the  spray. 
The  evaporation  rate  is  a  function  of  boiling 
point  and  surface  area  to  volume  ratio  (Asl 
V)  of  the  droplet.  The  evaporation  time 
increases  with  increasing  boiling  temperature 
and  SMD. 

2.  A  group  of  parameters  including  droplet  distri¬ 
bution,  evaporation,  mixture  formation  and 
subsequent  combustion  with  preheating  and 
dilution  of  reactants. 

3.  In  flameless  combustion  mode,  the  increased 
dilution  of  fresh  reactants  with  hot  combustion 
products  results  in  reduced  reaction  rate.  Due 
to  this,  reaction  zone  is  uniformly  distributed 
throughout  the  combustor  volume  with 
lower  peak  flame  temperature  than  that  of 
conventional  mode  [4].  In  conventional  mode, 
the  fuel  spray  directly  enters  the  combustion 
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zone  having  higher  peak  temperature.  There¬ 
fore,  the  droplet  evaporation  rate  is  relatively 
slower  in  flameless  combustion  mode  [4.  To 
achieve  complete  evaporation  and  combustion, 
the  droplet  residence  time  should  be  higher  in 
flameless  combustion  mode. 

To  sustain  flameless  combustion  with  liquid 
fuels,  the  above  three  issues  can  be  addressed  by 
increasing  the  residence  times  and  recirculation 
as  compared  to  flameless  combustion  with  gas¬ 
eous  fuel.  A  swirl  flow  creates  a  central  vortex 
zone  and  low  pressure  gradient,  supporting  a  large 
reverse  flow  region  in  the  combustor.  High  swirl 
creates  higher  centrifugal  force  that  enhances  the 
residence  time  of  the  hot  gases  trapped  within 
the  swirling  flow  [13-20].  The  increased  residence 
time  enhances  the  flame  stability  limits  and  rate 
of  mixing  of  products  and  reactants.  The  high 
recirculation  allows  for  good  mixing,  essential 
for  obtaining  a  distributed  reaction  zone  over  a 
large  volume  of  the  combustor  1 3—20  .  Therefore, 
a  tangential  air  injection  scheme  was  used  in  this 
study  to  generate  the  swirl  flow  in  the  combustor. 

In  this  study,  a  conical  combustor  with  60° 
diverging  angle  was  considered  with  a  total  vol- 
ume  of  ~0.004  m  [16,21  .  A  pressure  swirl  fuel 
injector  was  used  for  injection  of  kerosene.  The 
combustor  configuration  is  shown  in  Fig.  .  Com¬ 
putational  and  experimental  studies  were  carried 
out  simultaneously.  Recirculation  of  combustion 
products  was  identified  as  the  key  factor  to  sustain 
flameless  combustion.  Hence,  the  reactants  dilu¬ 
tion  ratio  (Rdii)  is  the  governing  metric.  Rdd  is  cal¬ 
culated  as  follows  [14,16  . 

p  ttlaxial  ( th0x  T  f ) 

( mox  +  mf) 

th  axial  py axial  dydz 

2.1.1.  Challenges  in  scaling 

Initially,  the  combustor  was  tested  at  21.5  kW 
thermal  input  using  the  unmodified  combustor 
configuration  (without  chamfer,  i.e.  Rc  =  0  mm) 
and  exit  diameter  (D)  of  25  mm.  Well  stabilized 
flameless  combustion  was  observed  experimen¬ 
tally  [16,21].  Computational  results  showed  that 
Rdu  varied  spatially  from  1.1  to  3.2.  The  same 
combustor  was  tested  at  higher  thermal  inputs 
of  40.8,  65.1  and  84.7  kW  (respective  heat  intensi¬ 
ties  10.2,  16.3  and  21.1  MW/m3)  and  experimental 
observations  revealed  that,  flameless  combustion 
was  not  stabilized  in  the  combustor  at  these  higher 
inputs  and  large  quantities  of  unburned  fuel  accu¬ 
mulated  in  the  combustor.  The  fuel  spray  cone 
angle  was  maintained  constant  at  45  for  all  noz¬ 
zles  and  SMD  of  all  nozzles  were  in  the  range  of 
35-37  pm  (details  in  Section  4.1).  Therefore,  with 
the  increased  fuel  mass  flow  rate,  droplet  number 


Fig.  1.  Dimensional  details  of  the  combustor. 


density  (DND)  also  increased.  Hence,  more  recir¬ 
culation  (increased  residence  time)  would  be 
required  to  increase  the  entrainment  and  to  achieve 
complete  evaporation.  The  computational  results 
also  revealed  that,  Rdd  decreased  with  increasing 
thermal  input.  Computational  and  experimental 
evidence  suggested  that  Rdd  needed  to  be  suffi¬ 
ciently  high  for  all  thermal  inputs  to  provide  the 
required  entrainment  and  residence  time. 

CRT  approach  appears  suitable  for  scaling 
swirl  combustors  operating  with  gaseous  fuels 
[12  .  However  in  case  of  liquid  fuels,  the  DND 
increases  with  thermal  input.  The  residence  time 
should  be  increased  for  complete  evaporation 
and  combustion.  CV  scaling  approach  for  higher 
thermal  inputs  results  in  increased  combustor  vol¬ 
ume  and  reduced  heat  intensity  [7,12].  Hence,  for 
the  present  case,  both  CV  and  CRT  approaches 
are  not  suitable.  Similarity  of  certain  dimension¬ 
less  quantities  in  scaling  is  different  for  various 
combustor  configurations,  operating  conditions 
and  modes  of  combustion  [22  .  Therefore,  a  com¬ 
bination  of  experimental  and  numerical  simula¬ 
tions  aimed  at  improving  the  droplet  residence 
times  and  recirculation  rates  were  considered.  To 
enhance  both  droplet  residence  time  and  recircula¬ 
tion  rate,  a  chamfer  at  the  top  of  the  combustor  is 
provided  as  shown  in  Fig  .  A  computational 
analysis  (described  below)  was  carried  out  for 
high  thermal  inputs,  by  varying  chamfer  radius 
( Rc )  to  determine  the  Run  for  each  case  and  the 
operating  conditions  that  yielded  an  Rdn>  2.5. 
The  combustor  geometry  for  different  thermal 
intensities  is  non-dimensionalized  with  D,  as 
shown  in  Fig.  1  and  listed  in  fable  2. 

2.1.2.  Numerical  method 

A  general  purpose  CFD  code  Fluent-14.5  was 
used  for  computational  studies  in  this  work.  A 
3-D  double-precision  pressure-based  solver  was 
used.  For  all  thermal  inputs,  tangential  air  inlet 
and  combustor  exit  velocities  were  maintained 
constant  to  ensure  similar  level  of  pressure  drop 
across  the  combustor.  Therefore,  the  air  inlet 
diameter  (din)  and  exit  diameter  ( D )  are  increased 
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Table  2 

Dimensional  details  of  the  combustor. 


Qth 

Q'" 

d 

D 

d\ 

h 

h 

Rc 

din 

21.5 

5.37 

5 

25 

3.2 

8.4 

4.4 

2.4 

0.4 

0.2 

40.8 

10.2 

7.5 

29 

2.758 

7.241 

3.793 

2.068 

0.689 

0.258 

65.1 

16.3 

9 

34 

2.325 

6.176 

3.235 

1.764 

0.735 

0.264 

84.7 

21.1 

11 

37 

2.162 

5.675 

2.972 

1.621 

0.81 

0.297 

with  increased  thermal  input.  Chamfer  radius  is 
varied  from  10  to  30  mm  at  5  mm  increments. 

Three-dimensional  Navier-Stokes  equations 
were  discretized  and  solved  in  a  finite-volume 
domain.  Reynolds  Stress  Model  (RSM)  was  used 
for  turbulence  modeling.  The  energy  equation 
was  solved  considering  20  intermediate  species 
equilibrium  chemistry  and  a  non-premixed  droplet 
combustion  model  for  simulating  the  combustion 
of  the  liquid  fuels.  Compressible  flow  was  consid¬ 
ered  and  the  viscosity  was  calculated  using  Suther¬ 
land’s  law.  Specific  heats  were  defined  as  a  function 
of  the  temperature  (piecewise-polynomial).  A  PI 
radiation  model  was  used.  Constant  mass-flow 
inlet  condition  normal  to  the  boundary  surface 
was  applied  at  air  inlets,  and  a  pressure  outlet  based 
boundary  condition  was  applied  at  the  exit.  No-slip 
wall  and  constant  temperature  boundary  condi¬ 
tions  were  applied  at  the  walls.  Non-premixed 
droplet  evaporation  and  combustion,  following 
the  spherical  law  was  considered  with  PDF  droplet 
evaporation.  A  single  component  surrogate, 
Ci2H23  was  used  to  simulate  kerosene  with  a  den¬ 
sity  of  780  kg/m3  Fuel  injection  was  simulated  as 
a  solid  cone  type  spray  with  a  droplet  diameter  of 
36  pm  and  a  cone  angle  of  45°.  The  amount  of  heat 
removal  from  the  combustor  walls  is  3.2,  8.3,  and 
12.9  kW  respectively,  for  three  higher  heat  intensi¬ 
ties  of  10.2,  16.3  and  21.1  MW/m3.  The  heat 
removal  through  wall  cooling  is  considered  by 
applying  heat-loss  through  combustor  walls  as  heat 
flux  boundary  condition  for  higher  heat  intensity 
cases.  The  solution  is  considered  to  be  converged 
when  RMS  residuals  of  the  system  were  less  than 
1  x  10  6.  A  number  of  computations  were  carried 
out  using  hexa  mesh  with  different  mesh  sizes  vary¬ 
ing  from  1.1  to  2.5  mm.  The  number  of  cells  for 
computations  was  varied  from  2  to  4.5  million 
elements.  A  mesh  size  of  1.2  mm  was  considered 
sufficient  to  obtain  grid-independent  results  with 
approximately  3.6  million  grid  points.  The  grid 
convergence  was  calculated  based  on  the  Grid 
Convergence  Index  ( GCI)  criteria.  If  the  GCI  for 
two  successive  grid  sizes  was  below  3%,  it  was 
considered  that  grid  convergence  has  been 
achieved  1  . 

2.2.  Reactants  dilution  ratio  ( Run) 

Run  was  calculated  at  different  axial  planes  of 
the  combustor.  For  the  case  of  g//7  =  21.5kW 


(Rc  =  0  mm)  with  exit  port  diameter  of  25  mm, 
a  maximum  Rdil  of  3.2  was  achieved.  Complete 
flameless  combustion  with  low  emissions  was 
observed  experimentally.  Rdd  was  calculated  for 
higher  thermal  inputs  (Qth  =  40.8,  65.1  and 
84.7  kW)  with  different  Rc  of  10,  15,  20,  25  and 
30  mm  and  results  are  shown  in  Ag.  2.  It  was 
observed  that,  for  constant  thermal  input  and 
increasing  Rc ,  the  degree  of  flow  reversal 
increased  in  the  combustor.  The  resulting  Rdd 
increased  with  Rc.  For  instance,  at  Qtjx  —  21.5- 
kW,  a  maximum  Rdd}  of  5.22,  5.77,  6.26  and 
6.75  was  obtained  for  Rc=10,  15,  20  and 
25  mm  respectively  (Fig.  2a).  The  curved  profile 
of  the  combustor  dome  and  chamfer  near  the  exit 
combined  to  form  curved  vanes  promoted  a  large 
degree  of  flow  reversal.  Hence,  Rddi  increased  with 
increased  Rc.  For  a  constant  Rc  and  increasing 
thermal  input,  Rdn,  was  calculated  as  shown  in 
Fig.  2a-d.  For  example,  at  Rc  =  25  mm,  the  max¬ 
imum  Recalculated  were  6.77,  3.75,  3.51  and  2.71 
respectively,  for  21.5,  40.8,  65.1  and  84.7  kW. 
Therefore,  it  was  observed  that  with  increasing 
thermal  input,  the  chamfer  radius  must  be 
increased  appropriately  to  maintain  a  constant 
Rdn  for  all  thermal  inputs.  The  zone  length  of 
Re>  2.71,  which  is  the  lower  limit  for  achieving 
flameless  combustion  is  calculated  for  all  compu¬ 
tational  conditions  and  shown  in  Ag.  2e.  It  was 
observed  experimentally  that  an  Rc  =  20,  25  and 
30  mm  were  sufficient  to  achieve  flameless  com¬ 
bustion  at  higher  heat  intensities  (10.2,  16.3  and 
21.1  MW/m3  respectively). 

2.3.  Residence  time  distribution 

The  residence  time  of  the  reactants  in  the  com¬ 
bustion  chamber  is  a  significant  parameter  to 
achieve  flameless  combustion  [23].  Three  basic 
time  parameters  were  considered  to  calculate  the 
residence  time  of  reactants. 

1.  Average  residence  time;Taug  =  (T);  V  —  com¬ 
bustor  volume  and  v=  volume  flow  rate  of 
reactants.  However,  since  the  present  combus¬ 
tor  operates  with  a  swirl  flow,  the  residence 
time  was  calculated  computationally  for  differ¬ 
ent  cases  by  injecting  particles  from  air/fuel 
inlets  in  the  combustor. 

2.  Swirl  based  residence  time  with  Rc  =  0  mm; 

^ s(Rc=0mm )  • 
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Fig.  2.  Variation  of  Rda  with  Rc  (a)  5.37  MW/m3  (b)  10.2  MW/m3  (c)  16.3  MW/m3  (d)  21.1  MW/m3  (e)  Minimum  and 
maximum  Rdit  for  all  cases  and  zone  length  of  Rdd>  2.72  (in  parenthesis). 


3.  Swirl  based  residence  time  with  Rc  =  25  mm; 

X s(Rc=25mm)  • 

The  calculated  xavg  decreased  from  0.61  to 
0.15  s  as  the  thermal  input  increased  from  5.37 
to  21.1  MW/m3.  The  calculated  t^c=o»™)  with 
swirl  flow,  decreased  from  0.76  to  0.29  s  and 
Ts(RC=25mm)  decreased  from  0.98  to  0.69  s  for  this 
same  range  of  thermal  inputs.  For  the  case  of 
21.1  MW/m3,  the  percentage  increase  in  residence 
time  is  93%  and  360%  for  Ts(Rc=0mm)  and  t s(Rc=25mm) 
respectively  as  compared  withiat;g.  It  was  observed 
from  the  computational  study  that  the  residence 
time  increased  with  both  swirl  flow  pattern  and 
increased  chamfer  radius  ( Rc ).  Residence  time 
distribution,  E(t)  [23]  was  calculated  for  all  cases 
by  considering  the  combustor  as  a  well-stirred 
reactor  and  shown  in  Fig.  .  The  E(t )  of  the  reac¬ 
tor  is  the  probability  density  function  of  a  particle 
in  the  reactor.  If  E(t)  of  a  reactor  is  high,  the  res¬ 
idence  time  of  the  particle  is  large.  It  is  observed 
from  Fig.  3  that  E(t)  increases  with  swirl  flow, 
and  increases  further  with  chamfer  plus  swirl  flow. 


3.  Details  of  experimental  methodology 

3.1.  Experimental  setup 

Figure  4  shows  a  schematic  diagram  of  the 
experimental  setup.  The  combustor  was  placed 
vertically  on  a  test  stand.  Kerosene  was  stored 
at  a  pressure  of  9  bar  (AP)  in  a  pressurized 
stainless-steel  tank.  The  fuel  injector  was  located 
at  the  center  of  the  combustor.  The  fuel  injector 
imparts  a  clockwise  rotation  to  fuel  spray;  hence 
a  counter-clockwise  air  injection  was  selected  to 
impart  more  shear  force  to  the  flow  resulting  in 
enhanced  mixing  and  evaporation  of  droplets. 
Air  supply  to  the  combustor  was  regulated 


-5.37  MW/m3  _xS(Rc  =  23mm)  =  0.98  s 
■21.1  MW/m3  Xsl^c  =  25  mml  =  0.69  s 
5.37  MW/m3  _tS(rc  =  o  mm)=  0.76  s 
21.1  MW/m3  _xs(Rc  =  0mm)  =  0.29  s 

5.37  MW/m3  _xAvg=  0.61  s 

10.2  MW/m3  xAvg=  0.29  s 

16.3  MW/m3  _xAvg=  0.20  s 
21.1  MW/m3  _  xAvg=  0.15  s 


Fig.  3.  Residence  time  distribution  (P(/))  for  three  cases 
based  on  x  s^Rc=Qmm^  and  xs^Rc=2s mm)  • 


Fig.  4.  Schematic  diagram  of  experimental  setup. 


through  electric  mass  flow  controllers  (accuracy 
±1.5%  of  full  scale). 

3.2.  Experimental  procedure  and  instruments 

Initially,  the  premixed  LPG-air  mixture  was 
ignited  with  a  spark  and  combustor  was  run  for 
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2-3  min  to  preheat  the  combustor.  The  kerosene 
fuel  is  injected  at  5  bar  pressure  by  opening  the 
ball-valve  in  the  fuel  line.  The  LPG  flow  rate 
was  then  gradually  reduced  and  the  kerosene 
injection  pressure  was  simultaneously  raised  to 
9  bar.  A  stable  flame  was  established  in  conven¬ 
tional  combustion  mode  with  stoichiometric  kero- 
sene-air  mixture  for  next  4-5  min.  After  an  initial 
start-up  time  of  7-8  min,  the  combustor  wall  tem¬ 
perature  reached  ~900  K.  A  chamfered  flange  was 
placed  at  the  top  to  effectively  reduce  the  exhaust 
port  diameter  from  90  mm  to  a  diameter  ( D )  for 
the  particular  heat  intensity  (  "able  2).  The  con¬ 
ventional  flame  then  gradually  shifted  to  a  flame¬ 
less  combustion  mode.  This  strategy  was  adapted 
to  understand  and  evaluate  the  effect  of  exit  port 
diameter  variation  on  transition  between  conven¬ 
tional  (90  mm)  and  flameless  combustion  mode 
(30  mm).  The  present  combustor  can  be  started 
with  the  top  components  in  place  for  a  real  prac¬ 
tical  application. 

Exhaust  gas  composition  was  measured  with  a  gas 
analyzer  which  included  02  analyzer  (0-25% 
range,  0.1%  accuracy),  CO  analyzer  (0-10,000  ppm, 
±5  ppm  accuracy),  NO  analyzer  (0-5000  ppm, 
±1  ppm  accuracy),  CxH y  analyzer  (0-50,000  ppm), 
and  C02  analyzer.  Temperature  measurements  were 
carried  out  with  R-type  (dJunction  =  1  mm)  thermo¬ 
couples.  The  sound  level  at  the  exit  (100  mm  away 
from  axis)  of  the  combustor  was  measured  for 
different  combustion  modes  with  a  fast  response 
(Resolution  =  0.1  dB,  t response  —  200  ms)  sound  level 
instrument. 


4.  Results  and  discussion 

4.1.  Spray  characteristics 

In  the  present  study,  four  nozzles  N1-N4  with 
mass  flow  rates  of  1.72,  3.27 ,  5.21  and  6.78  kg/h 
respectively,  were  used  to  provide  21.5,  40.8, 

65.1  and  84.7  kW  thermal  inputs  respectively. 
An  injection  pressure  of  9  bar  was  maintained 
for  all  experiments.  Various  details  of  the  spray 
characteristics  such  as  D10,  D32  (SMD),  Dvl0, 
Dv5o  and  Dy9o,  droplet  distribution,  droplet  num¬ 
ber  density  (DND)  were  measured  with  a  particle 
shadowgraphy  technique.  7000-9000  droplets 
were  considered  in  each  sample  size.  A  count  of 
150  pictures  was  selected  for  each  sample  at  an 
axial  position  of  45  mm  from  the  nozzle  tip. 

It  was  observed  that  for  all  four  nozzles,  SMD 
was  in  the  range  of  35-37  pm  and  variation  in 
other  diameters  was  relatively  very  small.  Since 
the  spray  cone  angle  and  droplet  diameters  were 
nearly  the  same  for  all  four  nozzles,  the  DND 
increased  for  higher  mass  flow  nozzles,  the  mea¬ 
sured  DND  for  N1-N4  nozzles  was  32  x  103, 
64  x  103,  110  x  103  and  167  x  103  n/cm3  respec¬ 
tively.  Therefore,  entrainment  of  hot  gases  needed 


to  be  increased  significantly  with  the  increasing 
DND  to  achieve  complete  evaporation  of  all 
droplets.  The  DND  distribution  for  all  four 
nozzles  is  shown  in  Fig.  . 

4.2.  Temperature  distribution 

Temperature  variation  in  the  radial  direction 
of  the  combustor  at  an  axial  location  of  120  mm 
was  measured  for  different  heat  intensities  at 
0  =  0.92  and  comparison  with  predicted  results 
is  shown  in  Fig.  6.  Due  to  larger  thermocouple 
response  time  (~0.25  s)  as  compared  to  integral 
turbulence  time-scales  (~3  ms),  it  is  difficult  to 
measure  actual  temperature  variation  in  the  com¬ 
bustor.  However,  temperature  variation  with  time 
is  measured  at  a  given  location  and  the  mean  was 
calculated  from  recorded  temperatures  over  a  per¬ 
iod  of  10-20  s.  The  measured  temperature  was 
corrected  by  considering  convection  and  radiation 
losses  from  the  thermocouple  junction.  For  the 
case  of  5.37  MW/m3,  the  wall  temperature  of  the 
combustor  was  ~800  K.  When  the  combustor 
was  operated  at  10.2  MW/m3,  the  walls  became 
red  hot.  Hence,  cooling  of  outer  walls  of  the  com¬ 
bustor  was  mandatory  for  higher  heat  densities, 
achieved  through  water  circulation  through  cop¬ 
per  tubes  brazed  on  the  outer  walls  of  the  combus¬ 
tor.  A  constant  wall  temperature  of  ~950  K  was 
maintained  for  higher  heat  intensities  (10.2- 

21.1  MW/m3).  The  heat  removal  through  wall 
cooling  is  3.2,  8.3,  and  12.9  kW  respectively  for 
three  higher  heat  intensities  of  10.2,  16.3  and 

21.1  MW/m3.  Fresh  air  at  ambient  temperature 
entered  the  combustor  and  circulated  on  the  inner 
walls;  a  sharp  rise  in  temperature  of  the  air  was 
observed  near  the  walls  of  the  combustor 
(Fig.  6).  Temperature  at  all  radial  locations 
increased  with  increasing  heat  intensity  of  the 
combustor.  As  expected,  the  temperature 
increased  from  the  walls  to  the  center  line  of  the 
combustor.  The  temperature  difference  across 
the  plane,  from  axis  to  near  wall  (0.0975  m)  for 
5.37  MW/m3  with  Rc  =  0  mm  was  443  K.  The 
temperature  difference  for  higher  heat  intensities 


Fig.  5.  Distribution  of  DND  for  nozzles  of  different  fuel 
flow  rates. 


V.M.  Reddy  et  al.  I  Proceedings  of  the  Combustion  Institute  35  (2015)  3581-3589 


3587 


Radial  location  (m) 


Fig.  6.  Temperature  distribution  comparison  of  exper¬ 
imental  and  computational  measurements  (-■-  5.37 
MW/m3_Rc  =  0  5.37  MW/m3_Rc  =  20  ••<>•• 

10.2  MW/m3_Rc  =  25  •  16.3  MW/m3_Rc  =  25 

21.1  MW/m3_Rc  =  30  ■  5.37  MW/ 

m3_Rc  =  20_Comp  — 10.2  MW/m3_Rc  =  25_ 
Comp  — ± —  16.3  MW/m3_Rc  =  25_Comp  — 

21.1  MW/m3_Rc  =  30_Comp. 

(10.2-21.1  MW/m3)  was  319,  293  and  245  K 
respectively.  With  increased  heat  intensity,  the 
overall  temperature  of  the  combustor  and  the 
temperature  of  the  fresh  air  circulating  increased. 
Hence  the  temperature  gradient  across  the  radial 
direction  decreased  significantly.  Maximum  tem¬ 
perature  at  the  center  of  the  combustor  increased 
from  1633  to  1741  K  as  heat  intensity  increased 
from  5.37  to  21.1  MW/m3.  The  temperature 
fluctuations  around  the  mean  value  were  in  the 
range  of  1.3-1. 8%  for  all  cases  (variation  bands 
shown  in  Fig.  6).  A  low  temperature  gradient 
and  smaller  fluctuations  are  representative 
characteristics  of  flameless  combustion.  For  all 
thermal  input  conditions,  the  maximum  tempera¬ 
ture  is  below  1800  K.  Therefore,  NO  v  emissions 
were  expected  to  be  relatively  very  low.  The 
predicted  temperatures  in  the  central  zone  are 
slightly  lower  than  the  measured  temperatures 
for  all  thermal  inputs.  For  the  outer  region  (next 
to  the  central  zone),  the  predicted  temperatures 
are  slightly  higher  than  measured  temperatures. 
Uniformly  distributed  temperature  with  low  tem¬ 
perature  gradients  is  observed  in  computational 
studies. 

4.3.  Pollutant  emissions 

The  CO,  NOx  and  HC  emissions  were  mea¬ 
sured  for  the  range  of  operating  conditions  and 
emission  levels  were  corrected  to  15%  02  level 
and  shown  in  Fig.  .  CO  emissions  increased  with 
a  decrease  in  0  from  1  to  0.6  and  increase  in  heat 
intensity.  However,  the  specific  emissions  index 
(ppm/kW)  decreased  with  increasing  heat  inten¬ 
sity.  For  Rc  =  20  mm  and  Q'"=5.31 ,  21.1  MW/ 
m3,  CO  emissions  varied  from  11  to  21  ppm  and 
25  to  41  ppm  respectively,  as  0  varied  from  1  to 


0.6.  The  specific  CO  emissions  for  these  cases 
varied  from  0.51  to  0.977  ppm/kW  and  0.3- 
0.48  ppm/kW  respectively.  The  emission  release 
rate  decreased  with  increasing  heat  intensity,  indi¬ 
cating  a  positive  outcome  for  higher  heat  density 
combustion  systems.  NOx  emissions  decreased 
with  decreasing  0 ,  as  expected.  For  lean  mixtures, 
the  average  measured  temperature  in  the  combus¬ 
tor  decreased  with  a  decrease  in  0.  This  led  to  a 
reduction  in  the  NOv  emissions,  however,  CO 
emissions  increase  slightly.  For  the  case  of 
Rc  =  20  mm  and  0"'=5.37.  21.1  MW/m3,  NO, 
varied  from  9  to  6  ppm  and  19-12  ppm  respec¬ 
tively,  for  0  varied  from  1  to  0.6.  The  specific 
NOx  emissions  for  these  cases  varied  from  0.42  to 
0.28  ppm/kW  and  0.22-0.14  ppm/kW  respectively. 

HC  emissions  increased  with  decreasing  0 
from  1  to  0.6  and  the  specific  emissions  decreased 
with  increasing  heat  intensity.  For  the  case  of 
Rc  =  20  mm  and  Q"= 5.37,  21.1  MW/m3,  HC 
emissions  varied  from  0  to  3  ppm  and  3-9  ppm 
respectively,  for  0  =  1-0.6.  The  specific  HC  emis¬ 
sions  for  these  cases  varied  from  0  to  0.14  ppm/ 
kW  and  0.03-0.1  ppm/kW  respectively.  The  over¬ 
all  variation  of  CO,  NOx  and  HC  emissions  for  all 
heat  intensities  (0  =  1-0.6)  were  measured  to  be 
ll^il,  6-19  and  0-9  ppm  respectively.  These 
emissions  are  well  within  the  range  of  emissions 
from  flameless  combustion  with  gaseous  fuels 
reported  in  the  literature. 

A  combustor  with  a  chamfer  radius  of 
Rc  =  25  mm  is  tested  for  all  thermal  inputs 
(21.5-84.7  kW)  conditions.  Flameless  combustion 
mode  is  observed  for  all  cases  without  any  issues 
related  to  combustion  stability.  Minimum  recircu¬ 
lation  required  for  each  case  of  thermal  input  is 
determined  experimentally  and  computationally 
by  varying  from  Rc  =  10-30  mm. 

4.4.  Acoustic  emissions 

Figure  8  shows  the  variation  of  acoustic  emis¬ 
sions  of  the  combustor  in  various  combustion 
modes.  Base  level  acoustic  emissions  of  84  dB 
were  measured  initially  for  cold  flow  conditions. 
After  ignition,  initially  the  combustor  operated 
in  the  conventional  mode  with  exit  diameter  of 
90  mm  and  the  level  of  acoustic  emissions 
increased  to  an  average  value  of  102  dB.  After 
3  min  of  conventional  combustion,  the  chamfered 
portion  was  mounted  and  the  exit  diameter  was 
reduced  to  D  mm  (  Table  2).  Immediately  after 
reducing  the  diameter,  the  sound  level  increased. 
After  a  time  of  2-3  min,  the  swirl  flow  was  well 
stabilized  in  the  combustor  and  flameless  combus¬ 
tion  was  observed.  The  sound  level  reduced  dra¬ 
matically  to  a  level  well  below  the  conventional 
combustion  mode.  For  the  case  of  21.1  MW/m3, 
113.5  and  93.6  dB  of  sound  level  was  observed 
in  the  transition  and  flameless  modes  respectively. 
It  was  observed  that  with  increased  heat  intensity, 
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Fig.  7.  Variation  of  emissions  with  equivalence  ratio  for  different  heat  intensities. 
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Fig.  8.  Variation  of  acoustic  emissions  for  all  heat 
intensities. 

the  sound  level  increased  during  the  operation  of 
the  combustor  in  transition  mode.  However,  for 
all  heat  intensities,  almost  a  same  sound  level  of 
approximately  94  dB  was  observed  during  the 
flameless  combustion  mode.  The  overall  net  sound 
level  reduction  from  conventional  to  flameless 
mode  for  all  combustors  was  in  the  range  of  8- 
9  dB.  A  similar  reduction  has  been  reported  in 
the  literature  7,14  . 
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5.  Conclusions 


profile  of  the  combustor  dome  and  chamfer 
combined  to  form  a  curved  vane  which  helps 
increase  the  degree  of  flow  reversal.  A  compu¬ 
tational  investigation  with  experimental  evi¬ 
dence  suggests  that  a  chamfer  radius  of  20, 
25  and  30  mm  was  sufficient  to  achieve  flame¬ 
less  combustion  for  g"'=10.2,  16.3  and 

21.1  MW/m3,  respectively. 

3.  The  peak  temperature  increases  in  the  combus¬ 
tor  and  the  temperature  gradients  decreases 
with  an  increase  in  the  heat  intensities.  The 
temperature  fluctuations  were  very  small 
(1.3-1. 8%  of  the  mean  value)  for  all  cases. 

4.  The  overall  variation  of  CO,  NOx  and  HC 
emissions  for  all  heat  intensities  (0  =  1-0.6) 
were  11^11,  6-19  and  0-9  ppm  respectively. 
These  emissions  are  well  within  the  range  of 
emissions  from  flameless  combustion  with  gas¬ 
eous  fuels  operating  at  high  intensity  in  the  lit¬ 
erature.  Specific  emissions  (ppm/kW)  decrease 
with  an  increase  in  heat  intensity. 

5.  The  outstanding  performance  of  the  burner 
with  very  low  chemical  and  acoustic  emissions 
at  high  heat  release  rates  indicate  the  potential 
for  use  in  various  industrial  and  gas  turbine 
applications. 


In  the  present  work,  a  new  combustor  configu¬ 
ration  was  designed  and  scaled-up  to  achieve 
flameless  combustion  with  liquid  fuels  at  high  heat 
intensities  for  various  industrial  and  gas  turbine 
applications.  Observations  are  summarized  below. 

1 .  Flameless  combustion  was  stabilized  in  the  base 
combustor  with  21.5  kW  thermal  input 
(5.37  MW/m3)  and  maximum  Rdm  of  3.2  with 
very  low  emissions.  However,  flameless  combus¬ 
tion  was  not  achieved  and  unburned  fuel  accumu¬ 
lated  in  the  combustor  for  higher  fuel  flow  rates. 

2.  A  chamfer  added  near  the  exit  in  the  modified 

combustor  configuration  helped  increase  the 
Ram  and  residence  time,  permitting  flameless 
combustion  at  higher  intensities.  The  curved 
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